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We perform a fully relativistic analysis of odd-type linear perturbations around a static and spher- 
ically symmetric solution in the most general scalar-tensor theory with second-order field equations, 
ft is shown that, as in the case of general relativity, the quadratic action for the perturbations re- 
duces to the one having only a single dynamical variable, from which concise formulas for no-ghost 
and no-gradient instability conditions are derived. Our result is applicable to all the theories of 
gravity with an extra scalar degree of freedom. We demonstrate how the generic formulas can be 
applied to some particular examples such as the Brans-Dicke theory, f(R) models, and Galileon 
,D ■ gravity. 
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I. INTRODUCTION 



General Relativity is definitely the most successful theory of gravity, passing all the experimental and observational 
tests performed so far. Nevertheless, it is still important to ask, "To what extent was Einstein correct?," and motivated 
by this question a number of alternative theories have been proposed. In particular, since the discovery of the cosmic 
acceleration [l[ , the mystery of dark energy has lead us to explore more seriously the possibility of modifying gravity 
on cosmological scales. 

In order for a given modified gravity model to be viable, it must pass stringent tests in the solar system and on the 
Earth. The standard cosmological evolution after nucleosynthesis must be reproduced as well in modified gravity. In 
addition to those requirements, one would demand that black hole solutions are stable, because there is strong evidence 
qq . that a black hole does exist e.g., at the center of our galaxy. The purpose of this paper is to provide the stability 
' criteria against linear perturbations about a static and spherically symmetric vacuum solution in modified gravity, 
, focusing, as a first step, on odd-parity modes. This is done by computing the quadratic action for the odd-parity 
f— ^ perturbations. We try to make our results as general as possible in such a way that the formulas we present are 
[ applicable to all the modified gravity theories having a single scalar degree of freedom cf> on top of the metric g^ v . To 
t— ( ■ this end, we employ the most general scalar-tensor theory with second-order field equations developed originally by 
Horndeski 0. Horndeski's theory contains novel models that have not ever been studied, as well as the Brans-Dicke 
• i-h . theory, f(R) models Q, and Galileon gravity 0, [1] as special cases. Therefore, our results are a useful and powerful 
discriminant to judge whether or not a given modified gravity model is viable according to the stability of a black 
H ! hole. 

We closely follow the approach of Refs. (HQ for black hole perturbations, in which two of the authors studied black 
hole perturbations in a gravitational theory whose Lagrangian is a general function f(R, C) of the Ricci scalar R and 
the Chern-Simons combination C . 

The paper is organized as follows. In the next section, we introduce the Lagrangian for the most general second-order 
scalar-tensor theory. After providing the background field equations in Sec. Ill, we derive in Sec. IV the quadratic 
action for odd-parity perturbations and a master equation for a single dynamical variable. From the quadratic action 
we give concise formulas for no-ghost and no-gradient instability conditions. Our generic formulas are then applied 
to some representative models of modified gravity in Sec. V. Finally, we draw our conclusion in Sec. VI. 
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II. THE MOST GENERAL SECOND-ORDER SCALAR-TENSOR THEORY 



We consider the most general scalar-tensor theory having second- order field equations both for the metric and the 
scalar field in four spacetime dimensions. The theory was originally determined by Horndeski and later it was 
rederived in Ref. Q as the most general extension of the Galileon field [3| • The equivalence of Horndeski's theory Q 
and the generalized Galileon Q was shown explicitly in Ref. . The generalized Galileon is described by the following 
four Lagrangians: 



C 2 = K{<j>,X), 

C 3 = -G 3 (0,X)D0, 

£ 4 = G 4 (cj),X)R + G iX - (V M V„0) 2 ] , 

£5 = G b {^X)G^V»V v ^-^G bX [{^f -ZU^V^V, 



2(V M V^) 3 ], 



(1) 

(2) 
(3) 

(4) 



where K and Gi arc arbitrary functions of 4> and X :~ — (d<fi) 2 /2. Here, we used the notation G%x for dGi/dX. Wc 
thus want to consider the action 



5 = E 



i=2 



d 4 x, 



(5) 



Horndeski's theory includes all the dark energy and modified gravity models with a single scalar degree of 
freedom in addition to the metric tensor. For example, C2 corresponds to the k-essence field. Since the non- 
minimal coupling of the form F{(f))R can be reproduced by choosing 6*4(0, X) = F((f>), the Brans-Dicke the- 
ory and f(R) gravity arc also special cases of the above theory. The covariant Galileon Q corresponds to 
K = -c 2 X,G 3 = -c 3 X/M 3 ,G 4 = M|j/2 - c 4 X 2 /M 6 , and G 5 = 3c 5 X 2 /M 9 , while the non-minimal coupling 
to the Gauss-Bonnet term, £(0) (i? 2 - 4R 2 ,, + R^ upa ), to K = 8£ (4) X 2 (3-lnX),G 3 = 4£< 3) X (7 - 3 lnX) , G 4 = 
4^X (2 - ]nX)jG 5 = -4fW lnX, where C (n) ■= d n ^/d<t> n @. The coupling of the form G^V ^ introduced 
e.g., in Refs. [To, |Tl| can be obtained from G5 = —(j) after integration by parts. 

Recently, Horndeski's theory has been revived and studied extensively in the context of cosmology [9l [l2| - [l6| . We 
are going to study Horndeski's theory from a different viewpoint, applying this powerful framework to the perturbation 
analysis in spherically symmetric spacetime. 



III. BACKGROUND 

Wc consider a static and spherically symmetric background, so that the unperturbed metric may be written as 

dr- 2 

g^dx"dx v = -A(r)dt 2 + — — — - + C(r)r 2 (d9 2 + sin 2 d<p 2 ) , (6) 

B(r) 

and the scalar field is dependent on the radial coordinate only: <f> = <p(r) (and therefore X = —B(<fr') 2 /2, where a 
prime stands for differentiation with respect to r.) Probably the easiest way to derive the background field equations 
is to substitute the above ansatz to the action and to carry out the variation with respect to A, B, C, and (j). Actually, 
the metric function C(r) is redundant, but the introduction of C(r) helps one derive the angular component of the 
gravitational field equations. Of course, we can set C(r) = 1 without loss of generality after getting all the field 
equations. The background field equations are thus given by 

£ A = 0, £ B = Q, £ c = 0, £ 4> = ^^^-{r 2 VABJ^-S = 0, (7) 

where explicit form of £a,£b, • • • etc. is provided in the appendix. The functions A(r), B(r) and <p(r) may be then 
determined by solving the field equations supplemented with appropriate boundary conditions. However, since we are 
interested in the generic structure of the perturbation equations and stability conditions rather than the properties of 
possible static and spherically symmetric solutions, we leave A, B, and <f> unspecified and proceed to the perturbation 
analysis. Note in particular that the background may not be an asymptotically flat solution. 
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IV. PERTURBATIONS 



Given the background equations of motion, we can now derive the quadratic action for perturbations. As is 
well known as the Regge- Wheeler formalism [Til , Il8j |. metric perturbations on a static and spherically symmetric 
background can be decomposed into odd- and even-type perturbations according to their transformation properties 
under a two-dimensional rotation. Such a decomposition is useful because the linear perturbation equations for odd 
and even modes are completely decoupled from each other, allowing one to study each mode separately. In this paper, 
we focus on the odd-type perturbations. We will obtain very concise formulas for the odd modes even in the most 
general case, while the computation of the even modes would be quite messy. A thorough analysis of the even-type 
perturbations will therefore be presented elsewhere [l9| . 

Each perturbation variable can be decomposed further into the sum of the spherical harmonics Yi m (6, ip). Then, 
at linear order in perturbations, or, equivalently, in the quadratic action, the perturbation variables having different 
£ and to do not mix, which simplifies our perturbation analysis drastically. Here we only summarize the definition 
and notation of the perturbation variables, referring the readers to Rcfs. p, @] for a more detailed explanation of the 
Regge- Wheeler formalism. 

The odd-type metric perturbations hn V = — can be written as 

h tt = 0, h tr = 0, h rr = 0, (8) 
hta = r)E ab d b Y im {6, if), (9) 

h ra = hi,tm(t, r)E ab d b Y em (8, tp), (10) 
h "b = \Y h Wrn{t, r) [E a c V c V b Y em {0 7 if) + E b c V c V a Y hn {6, <p)] , (11) 

where E ab := \J dct 7 e ab with 7^ being the two-dimensional metric on the sphere and e a b being the totally antisym- 
metric symbol with eo v = 1. The scalar <p do not have an odd- type perturbation. 

Because of general covariance, not all the metric perturbations are physical in the sense that some of them can 
be set to vanish by using the gauge transformation x 1 * — > .t' 1 + where £ M is an infinitesimal function. For the 
odd-type perturbations, one can completely fix the gauge by imposing the condition h% = for I > 2. This is called 
the Regge- Wheeler gauge, which we will use in the following calculations. For the dipole perturbation I = 1, h ab 
identically vanishes and we need to impose another gauge condition. We will therefore treat this case separately 
following the study for the general case. 



A. Second-order action 



Expanding the action (J5j to second order in the perturbations, we obtain the following action: 



5< 



2) 



J dtdrC {2 \ (12) 



where £' 2 ' is of the form 



^±lc {2) = ai + a 2 h\ + a 3 ( h\ - Ih^ + h' 2 + -hiho) , (13) 

Z7T \ r J 



where a dot denotes differentiation with respect to t. Since modes with different (£, m) do not mix, we have picked 
up the perturbations belonging to a given (£, m). Thanks to the spherical symmetry of background spacetime, the 
action for the m ^ modes takes exactly the same form as that for m = 0, and hence we may set m = without loss 
of generality Here and hereafter we abbreviate the subscripts £ and to. The coefficients 04,0,2, and 03 in the above 



4 



Lagrangian are given by 



"1 




(£a - £b) 



where 



T 
Q 
H 



2 ( G 4 + -Bct>'X'G bX - XG 5 ^ 



A' 

G 4 - 2XG iX + X ( — B<f>'G 5 x + G 54> 



G4 — 2XG4X + X \ G^x + G 



2A 88, 



c 



B OA" 



(14) 
(15) 
(16) 

(17) 
(18) 
(19) 



Note that the background field equations imply £a = = £b, allowing us to make these coefficients a little simpler. In 
the above equations, we have left £ a and £ b as they are in order to clarify the places where the background equations 
of motion are used. 

Notice that only G4 and G5 appear in the coefficients. This is also the case in tensor cosmological perturbations 
(corresponding to gravitational waves) at first and second order [H, [l4| . Given that the odd modes represent tensor 
degrees of freedom, it is reasonable to have such a similarity between perturbations on the different backgrounds. One 
may also notice that G5 appears only in the form of derivatives, which is consistent with the fact that £5 reduces to 
a total derivative if G5 = const. 

Since h is absent in the quadratic Lagrangian (|T5|) . ho is an auxiliary field. However, the presence of the term 
h[f hinders removing ho from the Lagrangian in a straightforward manner. To eliminate the non-dynamical degree of 
freedom, we invoke the following trick. Let us rewrite Eq. (|13|) as 



21+ 1 
2tt 



2a' 



,-a 3 



h + a 2 h\ + a 3 



-Q 



2q ( hi — h' + -ho 



(20) 



It is easy to verify that the new Lagrangian ((20)) reduces to the original one (fl~3)) after eliminating the auxiliary field 
q. The new Lagrangian enables us to put all the derivatives on q by integration by parts, transforming both ho and 
hi to removable auxiliary fields. Variation with respect to ho and hi now yields the following algebraic equations 
respectively: 



r [ra' 3 q + a 3 (rq' + 2q)] 



r 2 ai 



hi = —q. 

0-2 



(21) 

2(m 3 + a 3 ) a 2 

Plugging this result back into Eq. (|2"0"|) . we obtain the quadratic Lagrangian solely in terms of the master variable q. 
Once solving for q, one can determine the metric perturbations ho and hi through the relation (|21l) . 
To arrive at our final result, it is convenient to rescale the variable q as 



AT 



1/2 



(22) 



With this rescaling, the coefficient of the gradient term is transformed to —1 and the structure of the final Lagrangian 
will be of the clearest form. Here, we implicitly assumed that T > for the field redefinition (f2"2"j) to be feasible. 
Otherwise, the system would exhibit gradient instability. We thus obtain the quadratic Lagrangian for Q: 



2£+l 



£(2) 



1) 



2tt ~ 2(£-l)(£ 
where the effective potential V(r) is given by 

1 

"16 
-4 ( 3 



2) 



~ABQ 



Q 2 -Q 



12 



£{£ + \)T 
r 2 BU 



Q 2 - V(r)Q 7 



(23) 



V(r) 



4 




W 3A 


n'2 

+ 3 §r 






n + -rj 






2T" 


AT' 


2W 








rU + ~ 



B[_r_ 



7b 



~B 



10J" \ 
Br 2 H J 



(24) 
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Noting that we are imposing J- > in order to avoid gradient instability, we find a very concise expression for the 
no-ghost condition: 

Q > 0. (25) 
The equation of motion derived from the Lagrangian (|23[) is 

T ;Q - Q" + ^^IY q + V(r)Q = 0. (26) 



ABQ^ ^ r 2 BU 

Since t{t + 1) corresponds to the two-dimensional Laplacian in the real space, the third term in Eq. (]26f) represents 
wave propagation along the angular direction. Therefore, the squared propagation speeds along the radial direction, 
c^, and the angular direction, eg, are given respectively by 

c? = §, cg = |, (27) 

which in general differ from unity. The positivity of the first one is equivalent to the stability conditions mentioned 
above. To ensure the stability in the angular direction we further impose 

U > 0. (28) 

Interestingly, the speed of gravitational waves depends on the propagation direction in general second-order scalar- 
tensor theories. To summarize, T > 0, Q > 0, and H > are the necessary conditions for the stability of a static and 
spherically symmetric solution. 



B. Dipole perturbation: 1 = 1 

As we mentioned earlier, the discussion so far cannot be applied to the I = 1 mode. This may be also understood 
from the final Lagrangian (j2"3"]l . which is ill-defined for I = 1 due to the overall factor (t— For I = 1 case one 

can eliminate h\ by using the residual gauge degree of freedom. Furthermore, we note that ho contains a gauge mode, 
ho = C(t)r 2 , where C(t) is an arbitrary function of time. The equation of motion for ho is given by 

h'o - -ho = 0, (29) 
a s h'o' + a' 3 h' - ^^ho = 0, (30) 
with a3 = HyJ B/A. The solution is found to be 

where J and C(t) are integration constants. The second term is nothing but the gauge mode and hence can be elimi- 
nated by using the residual gauge degree of freedom. The first term cannot be removed by any gauge transformations 
and therefore is a physical perturbation. More specifically, it corresponds to a metric around a slowly rotating black 
hole and J is related to its angular momentum. This point is clearly understood by considering the case in general 
relativity: in that case the solution reduces to ho = — J/(47rM|>r), which coincides with the Kerr metric expanded to 
first order in the angular momentum J. 



V. APPLICATION TO PARTICULAR MODELS 



The results obtained in the previous section are so general that they can be applied to any scalar-tensor theory 
having second-order field equations. To demonstrate how useful our formalism is, in this section we will apply our 
stability conditions to some particular examples that have been frequently studied in the literature as potential dark 
energy and modified gravity models. 
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A. Non-minimally coupled models 

We begin with a familiar example of non-minimal coupling of the form F{<p)R. Since K and G3 are uncorrelated 
with stability against odd-parity perturbations, we may leave the two functions arbitrary. We are thus led to consider 
the following subclass of Horndeski's theory: 

K = K(X, 0), G 3 = G 3 (X, 0), G 4 = G 5 = 0. (32) 

This reduces to general relativity with a minimally coupled canonical scalar field for the choice K = X — V(<p), G3 = 0, 
and F = Mp^/2 = const, and to the original Brans-Dicke theory for the choice K = uiX/cf>, G3 = 0, and F = <fi with ui 
being a constant parameter. f(R) gravity can be mapped to a scalar-tensor theory and is reproduced with the choice 
K = f(cj)) - 0/'(0), G 3 = 0, and F = /'(</>). Kinetic gravity braiding d^-Hl corresponds to F{<j>) = with 
arbitrary K and G3. For this class of theories, we find 

T = Q=U = 2F(4>). (33) 

Therefore, no-ghost and no-gradient instability conditions are simply translated to 

F{<f) > 0. (34) 

The propagation speeds both along radial and angular directions are equal to the speed of light, c 2 . = c 2 , = 1. 



B. Covariant Galileon 

Covariant Galilcons correspond to 

X M 2 X 2 X 2 

K = -c 2 X, G 3 = - c 3^3, G 4 = - c 4 j^ , G 5 = 3c 5 -^g , (35) 

where C2, C3, c 4 , and C5 are dimensionlcss parameters and M is some mass scale. The stability conditions are given 
by 

F = M Pl + _(-<*+-— + -—) >0, (36) 

g = M n+M6{ c * + ^r) >0 > (37) 

9 6X 2 ( 2c 5 Bcb'\ , . 

n = M ^ + Jp{ Ci + "W J > °- (38) 

However, it was recentl y sh own that a static and spherically symmetric black hole cannot sustain a non-trivial profile 
of a covariant Galileon [2j|, implying that 0' = 0. Therefore, in the case of covariant Galileon the stability conditions 
are in fact trivial: Mpj > 0. 

C. Coupling to the Gauss-Bonnet term 

Let us consider a scalar field non-minimally coupled to the Gauss-Bonnet term, so that the gravitational part of 
the Lagrangian we want to study is 

M 2 

-^lR + ^4 > )(B 2 -ARl v + Rl vp(T ). (39) 

This Lagrangian is obtained from 

M 2 

Gi = -^+4^X(2-\nX), G 5 = -4£«lnX, (40) 
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where := d n ^/d(j) n , and we emphasize again that a concrete form of 7f and G3 is not important for stability 
against odd-parity perturbations. In this case the stability conditions read 

T = - 4£ (1) {B'4>' + 2B(f>") - 8^ (2) B0' 2 > 0, (41) 
Q = M^-A^^Bcj)' > 0, (42) 

H = M|, - 8£ (1) — > 0. (43) 
r 

As far as the authors know, the study of the odd-type perturbations around black-hole spacetime in this type of 
theories is provided only in Ref. [2~H . where a specific form for £(<f>) is assumed, £(4>) <x exp (<j>/ My{) . This form of 
£(4>) is, for instance, realized as the low-energy effective theory for the heterotic string. In Ref. [2J], the following 
master equation for the perturbations is derived: 

u"(r) + [V PC (r) +oj 2 K(r)] u(r) = 0, (44) 

where lo is a frequency of the wave. By analyzing this equation, the stability against perturbations is also established. 
However, K(r) given in Ref. p4| appears to depend on uj, which prevents us from checking consistency between our 
result and theirs. Let us also stress that the no-ghost condition is not provided in Ref. [24| and our stability conditions 
give further restriction on the scalar-Gauss-Bonnct coupling. 



D. Coupling to the Einstein tensor 

A scalar field derivatively coupled to the Einstein tensor has been studied recently in literature [ToL [TTI . |25| , for 
which the gravitational part of the Lagrangian is given by 

M 2 

-^-R + aG^d^d^, (45) 

where a is a parameter having dimension of (mass) -2 . This can be reproduced from 

M 2 

G 4 = ^i, G 5 = -a<f>, (46) 
after integration by parts. The stability conditions are given by 

T = Mp! - aBcj)' 2 > 0, g = H = Mp\ + aB<j)' 2 > 0. (47) 
The propagation speeds of gravitational waves along the radial and angular directions are therefore 



2 _ il/ 2 ! + aBcj)' 2 



A I 2 



4 = 1. (48) 



It can be seen that gravitational waves propagating along the radial direction are supcrluminal if a > and subluminal 
if a < 0. 



E. Possible application to higher-dimensional scenarios 

So far we have illustrated how one can apply our generic formulas to some representative models of modified gravity. 
Actually, the stability criteria we have derived can also be applied to rather different situations than this. Let us 
make a brief comment on a possible further application. 

Consider Lovelock gravity in D dimensions. The metric we want to discuss is of the form 

ds 2 = g^(x)dx»dx» + e 2 ^g mn (y)dy m dy n , (49) 

where g^ u (x) is the metric of four-dimensional static and spherically symmetric spacetime, and g m n(y) of (D — 
4)-dimensional compactified space. Since the Lovelock action yields second-order field equations for the metric, 
Kaluza-Klein reduction leads to a four-dimensional theory composed of <? M „ and <f> whose equations of motion are of 
second order. Therefore, the four-dimensional effective theory must be in a subclass of Horndcski's theory. This fact 
was emphasized in Ref. (26|. Our generic formula can therefore be used to investigate stability against odd-parity 
perturbations homogeneous along the extra dimensions. One of concrete examples is five-dimensional black strings 



in Gauss-Bonnet gravity; one can study the stability of the solution obtained in Ref. 27 | by using the results in this 
paper. 
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VI. CONCLUSION 



In this paper, we have studied odd-parity perturbations around static and spherically symmetric background space- 
time in the most general scalar-tensor theory with second-order field equations. Only tensor-type perturbations are 
excited in the odd-parity sector. It was shown, by explicit reduction of the quadratic action using the constraint 
equations, that there is a single master variable obeying a second-order differential equation, from which the dynam- 
ical behavior of all the perturbation variables are completely determined. From the quadratic action we derived the 
conditions to avoid ghost and gradient instabilities. They are summarized in a compact form as 

?(r) > °. S(r) > 0, H(r) > 0, (50) 

where the three functions are dependent upon the background metric and scalar-field profile, as defined in Eqs. (|17p 
(|19jl . Given a theory and a spherically symmetric configuration, one can judge by checking the above three conditions 
the stability of the spherically symmetric solution under consideration and whether or not the theory is viable. The 
squared propagation speeds in the radial and angular directions are given respectively by = QjT and <? e = G/T~L, 
which generically differ from the speed of light and even do not coincide with each other. This distinct feature in the 
propagation of gravitational waves is helpful to test some of alternative theories of gravity. 

In a subsequent paper [19| . we will report the analysis of linear perturbations in the even-parity sector. This can 
be done following the same line as in the present paper, but the computation is much more involved than in the 
odd-parity sector because the scalar field acquires an even-parity perturbation as well as the metric. It is reasonable 
to expect that the quadratic action for the even-parity sector contains two dynamical degrees of freedom, leading to 
two coupled second-order master equations. 
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Appendix A: Background equations 



Here we define the quantities that appear in the background equations |(7J). 



45 



£ A ■= K + B<j>'X'G 3X - 2XG 3<P + — {l-B- rB') G 4 + — [X'+- + —X G 4X + —XX'G 4XX 



X B' 



r B 



8B 



r A<j> 



AXG 



4<p<p ' 



2B0' ( X' X ) G i4>x 



[(1 - W)X' - 2B'X] G 5X 



' XX'G^ XX 



1 + B 



B')X + 2BX' 



4 , AB 

G50 Bcj) XG^sd, 

r r 



X' 



X 



XG 



5c/>X, 



(Al) 



£ B := K - 2XK X + + ^-J 'XG 3X + 2XG 3cl> 

+ H (iz£ - 4) Gt - i (izM _ 24) xc iX + ™ ( i 4) x>c iXX 

r \ r A J r \ r A J r \r A J 

4 A'\ /4 A'\ B6' A 1 
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£ c ■= K + 2XG H -X{B'ct>' + 2B<j)")G^x 



r V A \ V A / r 



G4 — 



B' 



2B(j)'' 



G. 



40 



X -5 



.A 



12 



A 2 



B' A'B 1 
2 — + —r- 

r A 



2B{A' + rA") 
rA 2 



Gax — ^^50 ] + BX' ( — \~ ~~ J (Gix — 



(rXA' + 2XA - rAX') G H x + 2XX' (2 + ^ ) 



2r 
BX 
r 



X 2 



A'B' 



A 
A' 



4' 2 A"B 
A 2 A 



2X' + — (rX 1 - 2X) 
A v ' 



G 



3A'B 
B 2 6'A' 



5<pX 



rA 



X' 

XX'G 5 xx 



r A 



4XG 



Atj>4> 



G5X - B<t>'X ( - + — ) 



(A3) 



J 



4 A 1 
r A 
A' 

XG^x 



XG 3X - 24>'G 3<t> + 20' 
1 - W A' 



1 - B B A' 



-2(j)' 



r A 
1-B B A' 



.4 

2B(j)' 



XG 



r 

2B A' 



G 



AX 



5A" 



X 2 G 



5XX 



S := 



r 2 r A 
—K^ + 2XG 3 ^ — Bcj)' X'G 3l px 



A[_ 
~A 



XG^x 



+B 



A[_ 
~A 



X' 



.A' 



A\ X 

~a)~ 



G 



4<t>X 



2A 



1-B B A' 

r 2 r A 



2B (1 
r \r 



A[_ 
~A 



XG 



AXX 



B/A'B' A A' A" 
~2 IT ~B + rl + ~A 



XG 



b<p<p 



X' + B—^ X ) G^^ 



(A4) 



G 



(A5) 
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